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Spirals  on  the  sea:  A manifestation  of  upper-ocean 
stirring 

Walter  Munk  and  Laurence  Armi 

Scripps  Institution  of  Oceanography,  La  Jolla,  California 


Abstract.  Spiral  eddies  were  first  seen  in  the  sun  glitter  on  the  Apollo 
Mission  30  years  ago;  they  have  since  been  recorded  on  SAR  missions  and 
in  the  infrared.  The  spirals  are  globally  distributed,  10-25  km  in  size  and 
overwhelmingly  cyclonic.  They  have  not  been  explained.  Under  light  winds 
favorable  to  visualization,  linear  surface  features  with  high  surfactant  density 
and  low  surface  roughness  are  of  common  occurrence.  We  have  proposed  that 
frontal  formations  concentrate  the  ambient  shear  and  prevailing  surfactants. 
Horizontal  shear  instabilities  ensue  when  the  shear  becomes  comparable  to  the 
Coriolis  frequency.  The  resulting  vortices  wind  the  linear  features  into  spirals. 
The  hypothesis  needs  to  be  tested  by  prolonged  measurements  and  surface 
truth.  Spiral  eddies  are  a manifestation  of  a sub-mesoscale  oceanography 
associated  with  upper  ocean  stirring;  dimensional  considerations  suggest  a 
horizontal  diffusivity  of  order  103m2s-1. 


Introduction 

The  first  photographs  of  spiral  eddies  appears  to 
have  been  taken  on  Apollo-Satum  in  October  1968.  In 
the  late  70s  SEASAT  with  its  synthetic  aperture  radar 
(SAR)  confirmed  the  early  discoveries  from  crewed 
space  flights  ( Stevenson  1998,  1999).  But  most  of  the 
existing  material  was  collected  by  Paul  Scully-Power 
(the  first  and  so  far  only  oceanographer-astronaut)  on 
5-13  October  1984:  “Far  and  away  the  most  impres- 
sive discovery...  is  that  of  the  submesoscale  ocean  (less 
than  100  km)  is  far  more  complex  dynamically  than 
ever  imagined...  Patterns  of  this  complexity  could  be 
seen  to  be  interconnected  for  hundreds  and  hundreds  of 
kilometers”  ( Scully-Power  1986). 

The  spiral  pattern  whose  global  distribution  was  re- 
ported by  Scully-Power  is  at  an  awkward  scale,  virtually 
impossible  to  recognize  from  shipboard,  and  too  large 
to  be  encompassed  even  from  high-flying  aircraft.  Dis- 
covery had  to  await  space  missions. 

Spiral  Images  in  the  Sun  Glitter  and  in 
SAR 

Figures  1 and  2 show  a visual  and  SAR  image,  re- 
spectively, of  spiral  streak  patterns.  Spirals  are  globally 
distributed  (Figure  3).  Typical  spiral  dimensions  are 
from  10  to  20  km,  with  streaks  50  to  100  m wide.  Spi- 
rals are  overwhelmingly  cyclonic,  wound  anti-clockwise 
(viewed  from  above)  in  the  northern  hemisphere,  clock- 
wise in  the  southern  hemisphere.  Ship  wakes  cross- 
ing the  streaks  (not  shown)  have  a cyclonic  offset  with 


Figure  1.  A pair  of  interconnected  spirals  in  the 
Mediterranean  Sea  south  of  Crete.  This  vortex  pair 
has  a clearly  visible  stagnation  point  between  the  two 
spirals,  the  cores  of  which  are  aligned  with  the  precon- 
ditioning wind  field.  7 October  1984. 

shears  up  to  10-3s-1.  We  need  to  refer  to  Munk  et  al. 
(2000)  (henceforth  MAFZ)  for  a more  representative  se- 
lection (13  images  out  of  400  collected). 

The  observational  material  poses  three  questions: 

(A)  How  are  the  spirals  wound? 

(B)  How  is  symmetry  broken  in  favour  of  cyclonic 
rotation? 

(C)  What  makes  spirals  visible? 

On  SAR  images  the  streaks  are  always  dark,  indicat- 
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Figure  2.  Spirals  in  the  Mediterranean  Sea  visual- 
ized with  Shuttle  XSAR.  The  streaks  are  differentially 
smooth.  9 October  1994. 
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Figure  3.  Distribution  of  spiral  eddies  from  Scully- 
Power  (1986)  visual  observations  and  our  collection  of 
400  images.  The  13  numbered  locations  refer  to  MAFZ. 
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Figure  4.  Rossby  number  in  the  upper  250  m sampled 
at  3 km  spacing  along  140°  W from  25°  N to  35°  N in 
the  North  Pacific  ( Rudnick  and  Ferrari  1999). 


ing  a reduced  scattering  cross-section,  e.g.  differentially 
smooth  water.  Natural  biogenic  surface  films  are  or- 
ganized by  near-surface  convergence  into  linear  streaks 
with  over  40%  surfactant  coverage  at  low  winds.  The 
concentration  is  associated  with  nearly  inextensible  sur- 
face films  which  dissipate  capillaries  and  short  gravity 
waves.  The  film  thickness  required  to  dampen  the  short 
waves  is  only  0.01  to  0.1  /xm.  On  the  optical  images  the 
smooth  streaks  are  bright  in  the  inner  sun  glitter  (which 
requires  low  rms  slopes  for  reflection  of  the  sun  into  the 
camera)  and  dark  in  the  outer  glitter.  The  situation 
is  complex  and  not  well  understood,  and  we  refer  to 
MAFZ  (1225-30,  1236-7)  for  further  discussion. 

But  evidently  the  third  question  can  be  restated  as 
follows:  What  is  the  circulation  pattern  that  collects  the 
surfactant  material  into  streaks  (which  are  subsequently 
wound  into  a spiral  pattern)?  Multiple  stripes  at  kilo- 
meter spacing  presumably  are  associated  with  helical 
circulation  rolls  in  the  atmospheric  boundary  layer.  In 
addition,  frontal  instabilities  can  concentrate  and  dis- 
tort the  surfactant,  as  we  shall  see. 

Ambient  Ocean  Vorticity 

Measurement  of  surface  velocity  shear  du/dy  along 
1000  km  of  roughly  northward  track  in  the  North  Pa- 
cific ( Rudnick  and  Ferrari  1999)  indicate  values  of  order 
10~5s_1.  The  situation  is  conveniently  portrayed  by  a 
distribution  of  Rossby  Numbers 

R°  = c//  (!) 

where  £ = dv/dx  — du/dy  is  the  vertical  component  of 
vorticity  (cyclonic  is  positive)  and  / is  the  Coriolis  fre- 
quency (Figure  4).  The  distribution  is  symmetric,  with 
very  few  values  exceeding  |.  There  are  a few  outliers 
showing  a slight  preference  of  cyclonic  vorticity  for  large 
|Ro|,  and  this  has  since  been  confirmed  ( Rudnick , this 
volume) . 

The  above  paper  also  shows  that  the  shear  is  dis- 
tributed over  a broad  band  of  scales,  from  kilometers  to 
hundreds  of  kilometers. 

Horizontal  Shear  Instability 

Starting  from  parallel  shear  flow  with  an  inflection 
point,  Figure  5 shows  a numerical  simulation  of  the 
development  of  the  most  unstable  mode  ( Corcos  and 
Sherman  1976,  1984).  Time  is  in  units  of  the  initial 
reciprocal  shear  at  the  stagnation  point  The  numerical 
experiment  was  intended  to  model  a vertical  shear  flow, 
but  may  as  well  be  interpreted  in  terms  of  a horizontal 
shear  flow.  There  is  no  implication  of  the  sense  of  ro- 
tation; in  fact  we  have  reversed  the  published  drawing 
from  anticyclonic  to  cyclonic  rotation. 
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Figure  5.  Computer  simulations  of  a developing  shear 
instability  ( Corcos  and  Sherman  1984).  The  four  pan- 
els show  the  streamlines  at  times  0.5,  1.0,  1.5,  2.0  (in 
units  of  the  initial  reciprocal  shear).  Heavy  line  is  the 
“cat’s-eye”  streamline  through  the  stagnation  points. 
The  dots  represent  particle  positions  initially  placed  on 
the  interface;  they  are  initially  crowded  near  the  two 
stagnation  points  to  allow  for  a subsequent  large  strain. 
The  model  allows  for  diffusion  and  viscosity.  We  have 
reversed  the  original  figure  from  anticyclonic  to  cyclonic 
rotation. 


readily  shown  that  vorticity  is  conserved: 

DC/Dt  = 0,  C = -V2tf  (3) 

where  DjDt  is  the  substantial  derivative  (following  a 
particle).  The  circulation  T = J ( dx  dy  = 4n  U k~x  is 
the  same  for  all  four  panels.  The  streamlines  through 
the  two  stagnation  points  separate  the  “core  circula- 
tion” from  an  exterior  circulation  which  is  not  all  that 
different  from  the  case  a = 0.  The  fractional  core  cir- 
culation 

Tc/ T = 47r~1tan~1  -y/a  (4) 

increases  with  increasing  a.  The  four  panels  give  in- 
dependent steady-state  solutions  for  different  values  of 
the  independent  parameter  a,  but  we  are  tempted,  by 
comparison  of  the  two  figures,  to  interpret  the  Stuart  so- 
lution as  a developing  instability  for  a slowly  increasing 
a(t).  In  MAFZ  we  demonstrate  that  an  initial  particle 
injection  produces  mature  spirals,  as  in  the  Corcos  and 
Sherman  numerical  experiment,  some  of  them  resem- 
bling the  vortex  pair  in  the  Cretan  Sea  (Figure  1).  The 
proposed  response  to  question  A is  that  the  spirals  are 
a manifestation  of  a horizontal  shear  instability. 

There  is  as  yet  no  breaking  of  symmetry.  In  fact, 
replacing  D^/Dt  = 0 by 

D(C  + f)/Dt~0  (5) 

does  not  change  the  developing  streamline  pattern,  so 
that  planetary  vorticity  is  conserved  in  the  /-plane. 
But  the  developing  pressure  patterns  are  fundamentally 
modified  ( MAZD  figure  22).  Starting  with  the  Stuart 
solution  of  a low  pressure  in  the  core  (regardless  of  the 
sense  of  rotation)  and  gradually  turning  on  the  corio- 
lis  parameter  /,  the  core  pressure  remains  a LOW  for 
the  cyclonic  case  but  switches  to  a HIGH  for  the  an- 
ticyclonic case.  Clearly  there  is  here  the  dynamics  for 
breaking  symmetry. 

Breaking  Symmetry:  the 
Hoskins-Bretherton  Jet 


Streamlines  show  the  development  of  Kelvin’s  cele- 
brated “cat’s-eye”  solution.  Particles  inserted  along  the 
interface  exhibit  the  growth  of  a spiral. 

It  is  instructive  to  compare  this  numerical  experi- 
ment with  Stuart’s  (1967)  steady-state  solution  for  a 
non-rotating  incompressible  equal-density  fluid  (Figure 
6).  The  streamlines  are  given  by 

'S?(x,y)  = — &-1f7  logfcosh  ky  — a cos  kx).  (2) 

For  a = 0 we  have  u = —d^f/dy  = -U tanh  ky.  It  is 


Hoskins  and  Bretherton  (1972)  have  solved  a prob- 
lem of  frontogenesis  with  conservation  of  density  and 
potential  vorticity, 

D{p,q)/Dt  = 0,  pg  = (/  + C)*Vp  (6) 

where  D/Dt  is  the  substantial  derivative.  The  start- 
ing point  is  a vertically  mixed  layer  with  a horizontal 
density  transition  from  warm  and  light  in  the  south 
(say)  to  cold  and  heavy  in  the  north  (Figure  7,  left). 
The  initial  density  gradient  develops  into  an  eastward 
“thermal  wind”,  as  shown.  A deformation  field  7 = 
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a = 0 a = 0.25 
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Figure  6.  Stuart’s  (1967)  solution  portraying  a steady 
flow  conserving  vorticity.  Streamlines  designate  both 
dimensional  streamlines  4'  = — logg,  q = cosh  y — 
acosx),  and  lines  of  constant  vorticity  ( = (1  —at2)q~2, 
with  values  indicated  in  each  panel  to  the  left  and 
right,  respectively.  For  a = 0 the  velocity  profile  is 
u = — tanhy.  Stagnation  points  at  x = ±ir,  y = 0 
are  connected  by  the  stagnation  stream  lines  separat- 
ing the  closed  “core”  circulation  from  the  outer  flow. 
With  increasing  a a larger  fraction  of  the  circulation  is 
concentrated  within  the  core. 

du/dx  — dv/dy  is  superposed,  causing  the  initially  ver- 
tical isopycnals  to  tilt  northward.  Up  to  this  point  there 
has  been  no  breaking  of  symmetry,  all  directions  can  be 
reversed. 

In  the  subsequent  development  we  use  “north”  and 
“east”  only  for  convenient  reference  to  the  figure.  The 
northward  tilting  of  isopycnals  is  not  uniform;  the 
northern  isopycnals  converge  at  the  surface,  and  the 
southern  isopycnals  diverge.  Accordingly  the  associated 
eastward  thermal  jet  has  a strong  cyclonic  shear  at  its 
northern  (left)  flank  and  a weak  anticyclonic  shear  at  its 
southern  (right)  flank.  At  time  2.5  (measured  in  7-1 
units)  the  associated  Rossby  numbers  are  Ro"r  = +1 
and  Ro~  = —0.3,  respectively. 

At  this  stage  the  underlying  rate  of  strain  no  longer 
determines  the  rate  of  development.  Rather,  the  isopy- 
cnal  ’’collapse”  takes  the  form  of  a “Rossby  Adjust- 
ment Problem”  with  /-1  « 10_4s_1  taking  the  place 
of  7-1  « 10-5s_1  as  the  relevant  time  scale  ( Ou  1984). 
In  the  short  time  interval  between  2.5  to  2.75  7-1  the 
cyclonic  shear  grows  from  Ro+  = 1 to  Ro+  = 3,  and 


at  2.89  7-1  the  density  gradient  at  the  left  flank  de- 
velops Ro+  = 00,  while  anticyclonic  shear  remains  at 
Ro-  = —0.3. 

The  crucial  point  is  that  starting  at  a time  when  Ro+ 
is  of  order  +1  the  cyclonic  shear  zone  becomes  a breed- 
ing ground  for  spiral  eddies  long  before  appreciable  an- 
ticyclonic vorticity  has  been  generated.  In  Figure  4 the 
third  panel  has  been  emphasized  because  at  this  time 
the  vertical  shear  at  the  surface  at  the  front  reaches  a 
value  of  du/dz  = 2N  (N  is  the  buoyancy  frequency)  cor- 
responding to  a Richardson  number  (N/ (du/dz))2  — ~ 
and  suggesting  the  onset  of  vertical  shear  instability. 

An  independent  consideration  has  to  do  with  the  vis- 
ibility of  the  spiral  arms,  presumably  the  result  of  the 
alignment  and  concentration  of  surfactants.  Consider 
an  elementary  surface  area  SxSy  at  time  zero.  With  the 
developing  front  the  area  is  elongated  along  the  x-axis 
on  both  flanks  of  the  developing  jet.  But  in  accordance 
with  the  Hoskins  and  Bretherton  theory,  at  the  time 
2.75  7-1  the  area  has  expanded  (by  a factor  |)  on  the 
anticyclonic  side,  while  it  has  contracted  (to  i the  orig- 
inal area)  on  the  cyclonic  side.  Thus  the  frontal  theory 
has  the  elements  to  account  for  both  the  visibility  and 
sense  of  rotation  of  the  spiral  eddies.  But  when  exam- 
ined in  detail  the  story  is  not  as  clear-cut  as  presented 
here,  and  we  must  refer  to  MAFZ  for  a more  complete 
discussion. 

Breaking  Symmetry:  the  Margules 
Front 

In  the  traditional  Margules  treatment  of  a sharp 
front,  anticyclonic  fronts  are  ruled  out  by  the  condi- 
tion of  static  stability.  For  a “softened  Margules  front” 
(MAFZ)  there  is  a gravitational  stability  limit  to  the 
anticyclonic  shear,  but  no  such  limit  to  the  cyclonic 
case. 

Breaking  Symmetry:  Anticyclonic 
Instabilities 

To  further  confuse  the  issue,  there  is  another  inde- 
pendent set  of  processes  to  explain  the  dominance  of 
cyclonic  vortices.  It  follows  from  the  Rayleigh  crite- 
rion of  stability  (extended  to  include  Coriolos  accel- 
eration) that  cyclonic  circular  vortices  are  stable  and 
anticyclonic  vortices  are  unstable  (MAFZ  6.21),  and 
this  leads  to  an  “inertial  instability”  criterion  Ro  < — 1 
which  goes  back  to  Pedley  (1969).  Oceanographers  are 
familiar  with  the  vertical  shear  (Richardson)  instability 
for  du/dz  > 2 N but  surprisingly  unfamiliar  with  the 
horizontal  shear  (Pedley)  instability  for  du/dy  > f. 

There  is  ample  numerical  evidence  for  instabilities 
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Figure  7.  Cartoon  for  the  generation  of  ocean  spirals 
(see  text). 


that  impede  the  development  of  anticyclonic  vortices 
but  have  no  such  effect  on  cyclonic  vortices.  ( Lesieur  et 
al,  1991,  Potylitsin  and  Peltier  1998).  The  distinction 
is  equally  clear  in  laboratory  experiments  ( Bidokhti  and 
Tritton  1992).  Kloosterziel  (1990)  observed  that  “it  is 
virtually  impossible  to  create  anticyclonic  vortices  ...by 
simply  stirring  the  fluid  locally.  But  cyclonic  vortices 
are  easily  created  in  this  way!  In  the  anticyclonic  case 
stirring  ....  leads  to  turbulent  motion  and  the  genera- 
tion of  waves  whereas  by  cyclonic  stirring  a well-defined 
smooth  vortex  forms.” 

With  regard  to  question  B,  it  is  unresolved  whether 
the  dominance  of  cyclonic  vortices  is  associated  with 
a dominance  in  cyclonic  horizontal  shear  early  in  the 
formation  process,  or  with  the  relative  instability  of  an- 
ticyclonic vortices  in  the  mature  stage. 


Stirring  and  Mixing 

The  proposed  formation  of  spiral  eddies  is  cartooned 
in  Figure  7.  The  geostrophically  balanced  ambient 
ocean  vorticity  of  order  ±10 _1/  is  enhanced  by  local 
frontogenic  processes,  with  concentration  of  surfactant 
along  a converging  line.  When  the  frontal  shear  be- 
comes comparable  to  f,  instabilities  lead  to  cross- frontal 
flow  accompanied  by  a cat’s-eye  circulation  pattern. 
The  cat’s-eye  circulation  twists  the  convergence  line  and 
neighboring  linear  features  into  a cyclonic  spiral  which 
stretches  and  further  thins  the  lines  of  surfactant  con- 
centration. It  is  a text-book  example  of  horizontal  stir- 
ring. 

We  have  not  made  any  serious  attempt  at  estimat- 
ing the  relevant  diffusivity.  The  best  we  can  do  is  a 
dimensional  argument  for  k = cl' u'  with  c = \ (typical 
for  turbulent  shear  flows,  Tennekes  and  Lumley  1972). 
Taking  the  minor  axis  of  a Stuart  core  as  the  charac- 
teristic length  V,  the  maximum  velocity  along  the  stag- 
nation streamline  for  u'  and  a — 1 for  a mature  20  km 
spiral  yields  k — 103m2  s-1  ( MAFZ ). 


Testing  the  Hypothesis 

Our  hypothesis  is  based  on  observational  material 
which  consists  almost  entirely  of  unrelated  glimpses  in 
x,  y- space  on  the  sea  surface.  For  a satellite  in  a low 
earth  orbit  (LEO)  a given  point  remains  within  view 
for  only  about  6 s.  What  is  required  here  are  prolonged 
stares  or  frequent  repeat  visits  coordinated  with  ship- 
board observations  (Figure  8).  We  cannot  think  of  any 
x,y,z,t  ocean  processes  that  had  been  properly  iden- 
tified from  measurements  in  half  the  coordinate  space. 
We  must  assume  that  there  are  serious  flaws  in  the  fore- 
going presentation. 


Figure  8.  The  proposed  experiment.  SAR  imagery 
from  an  overhead  drone  is  examined  by  the  authors  on 
shipboard  in  real  time.  The  image  shows  the  position 
of  the  vessel  which  is  about  to  enter  a spiral  streak. 

Following  his  41-G  space  mission  in  October  1984. 
Scully-Power  (1986)  wrote:  “The  almost  ubiquitous  oc- 
currence (of  spiral  eddies),  whenever  submesoscale  dy- 
namics was  revealed  in  the  sun  glitter,  indicates  that 
they  are  perhaps  the  most  fundamental  entity  in  ocean 
dynamics  at  this  scale.  The  difficulty  is  in  explaining 
their  structure.”  The  only  serious  attempt  at  analysis 
has  been  in  a Norwegian  Doctoral  dissertation  which  ex- 
plores baroclinic  instabilities  in  a narrow  cyclonic  shear 
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zone  ( Eldevik  and  Dysthe  1999).  Why  has  the  problem 
received  so  little  attention  in  the  thirty  years  since  dis- 
covery? We  assert  that  the  fashion  during  these  years 
has  been  statistical  rather  than  phenomenological  de- 
scriptions of  ocean  features,  and  here  we  are  concerned 
with  a truly  phenomenological  problem. 
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